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We investigate in detail recent suggestions that codimension-two braneworlds in 
six dimensional supergravity might circumvent Weinberg's no-go theorem for self- 
tuning of the cosmological constant. The branes are given finite thickness in order to 
regularize mild singularities in their vicinity, and we allow them to have an arbitrary 
equation of state. We study perturbatively the time evolution of the solutions by 
solving the equations of motion linearized around a static background. Even allowing 
for the most general possibility of warping and nonconical singularities, the geometry 
does not relax to a static solution when the brane stress-energies are perturbed. 
Rather, both the internal and external geometries become time-dependent, and the 
system does not exhibit any self-tuning behavior. 



INTRODUCTION 



There has been considerable recent interest in the possibility that braneworld construe 



m 



tions give new possibilities for solving the cosmological constant problem 
y,y, 0,0,^^. One rationale is that intrinsically extra-dimensional effects might provide a 
loophole to Weinberg's no-go theorem against self-tuning mechanisms, which is formulated 
in 4D [l^ . Such attempts first arose using 5D models, where solutions could be constructed 
with the propertythat, regardless of the 4D cosmological constant (brane tension), the uni- 
verse was static [lllSjIsl. Unfortunately it was shown that these solutions had singularities 
in the bulk which when regularized corresponded to the presence of an additional fine-tuned 
brane; moreover these solutions had flat directions which would give rise to nonstatic solu- 
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tions when the branes were perturbed away from their fine-tuned values. Other attempts 



15|, 
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involving black holes in the bulk were also plagued by instabilities ^ , 

The_case of a_6D bulk with (codimension-two) 3-branes has proved to be somewhat more 



id. 




11, 12, 



2L 



J2l|j2aj2 



2a, 



24 



Q, 22, 3,0,0,0,0, 



subtle 

E , 3, la^, llll, 1221, lia |41|, 142, 14^ |4^^ The conical singularity induced in 

the 6D Ricci scalar by the codimension-two source exactly cancels the contribution of the 
delta-function source itself when the 6D action is dimensionally reduced. The contribution 
of the brane tension to the 4D cosmological constant thus disappears automatically; the 4D 
effective action does not see it — it is as though it has been fine-tuned out of the theory, but 
in fact it was the dynamics of 6D general relativity which canceled it. Superficially, this 
looks exactly like the kind of inherently extra-dimensional effect which might well circumvent 
Weinberg's theorem. In practice however, the idea has not yet been shown to work. The 
problem is that the tension through its gravitational couplings affects other terms in the 
effective action besides the bare delta function source. In the models that have been studied 
so far, the sensitivity of these extra terms to the brane tension spoils the hoped-for self- 
tuning effect that would manifest itself when the tension undergoes a dynamical change, as 
during a phase transition. 

The failure of self-tuning in the simplest model, where the extra dimensions are compact- 
ified by the fiux of a U(l) gauge field, was demonstrated in references 12 1 and 41 1. In |4ll |. 
it was also argued that the same observations rule out self-tuning in supergravity versions of 
the model. We disagree with this point of view, as explained below. Confiicting arguments 
in the literature hold out the hope that self-tuning could still work in 6D SUGRA models 



44l |. so it is worthwhile to investigate in detail, given the notorious intractability of the 



cosmological constant problem. 

We review the formalism of codimension-two braneworld models in section II, with and 
without supersymmetry, and the details of arguments for and against self-tuning which have 
been made for the SUSY case. We also explain why it is useful to consider finite-thickness 
branes. In section III we set up and solve the leading-order equations of motion for the SUSY 
case, treating the time-dependent brane sources as perturbations. Section IV focuses on a 
particular class of warped background solutions which is crucial to settling the self-tuning 
issue; we demonstrate that these solutions are adequately described by our perturbative 
treatment so that self-tuning, if present, would be visible within our calculational framework. 
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In Section V we impose the boundary conditions at the branes, which reduces the general 
solution of section 111 to the particular ones of cosmological interest, in the case of time- 
dependent change of brane tension. We show that the system does not relax back to a static 
solution, but rather evolves like a universe with nonvanishing cosmological constant — there 
is no self-tuning. These results are discussed in section VI, and technical details are given 
in the appendix. 

II. REVIEW OF CODIMENSION-TWO BRANEWORLDS 

In this section we will briefly review how self-tuning is known to fail in the nonsupersym- 
metric 6D models, and the issues surrounding whether it can work in SUGRA versions of 
the models. 



A. Einstein-Hilbert gravity 



The simplest model which allows for compactification of the two extra dimensions has 
the action 



S 



~2 



1 

A' 



+ ^1 



branes ; 



where 



branes 



(1) 



(2) 



The interplay between a positive bulk cosmological constant and a gauge field flux Fro gives 
an energetically preferred size for the compact space. (It should be noted that this type of 
model was initially studied in the mid 1980 's, before the concept of branes was introduced 

Solutions exist with static extra dimensions and which are maximally symmetric in the 
large four dimensions, [53 1 



ds^ = -de + e^'''dx' + k-''(dr'' + sm^rd9^) 



A, 



^ n + ^ 
— (1 ± cosr) 



(3) 



with /c^ = M^*{Kf^/2 + 3/5V4), = M^^{Ae/6 - the coordinates r and 9 range 
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over (0, it) and (0, 27r) respectively, and where the 6D Ricci scalar has the form 

Thus, the singular parts of the action cancel exactly when evaluated for the classical so- 
lutions, and the effective 4D cosmological constant is determined entirely by the smooth 
parts: 

A4 = V, (1a6 - , (5) 
where V2 is the volume of the compact manifold, 

V2 = 2 (27r - 6) k'^, (6) 

and 6 is the deficit angle associated with conical singularities induced by 3-branes, so that 
in the presence of brane tension, the coordinate 6 defined above ranges from to 27r — 6. 
It requires a fine tuning of the bulk quantities Ag and (3 to have a vanishing cosmological 
constant, so the cancellation of the singular parts would not in itself explain the smallness 
of A4. Nevertheless it is suggestive; for example if supersymmetry in the bulk caused (0) to 
vanish, then SUSY could be broken on the branes without affecting A4. 

However this seemingly elegant idea does not ^york, as was explicitly shown in The 
reason can be understood intuitively as follows [41]. The simplest version of the model 
has two branes with equal tension r at antipodal points of the spherical extra dimensions, 
which produces conical defects and a deficit angle 6 = t/Mq. The volume of the extra 
dimensions depends linearly on the deficit angle, eq. ©• The flux of the gauge field in the 
extra dimensions is proportional to the volume, $ = Fr0V2- This flux is conserved, so if one 
tries to change the brane tensions, the field strength has to change to compensate for the 
change in volume. Any such change in Fro spoils the tuning needed to make A4 vanish in 
((HI). Thus A4 is linearly dependent on the brane tensions — not through the delta function 
terms in the original action, but as a consequence of the conical defect which they produce 
on the bulk geometry. 
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B. Super gravity 

A supergravity version of the above model involving chiral fields $a and a dilaton can 
be derived from the action 



S 



-e^t;(<l>) 



•S'branes) (7) 



One can consistently assume that the potential is minimized when the fields = 0, 
and thus ignore these fields and treat f ($) as the bulk cosmological constant, Ag. 

Comparing with the nonSUSY model, it is easy to see that there exist solutions with 
constant (/), where everywhere we have replaced Ag Age*^ and /^^e"*^, provided that 

takes the value 0o such that its induced potential is stationary, 

^'(0o) = -^/9'e-*'+Age*« = O (8) 

Comparing to (0), this is precisely the value needed to make the effective cosmological 
constant vanish, 

A4 - V2 (^Ag e-^" - -^(3^ e--^") = (9) 

Jiis seemingly automatic adjustment of A4 to zero is one reason that previous authors 
T, 36, 42 , Q, Q have been encouraged to seek a self-tuning mechanism in this model. 
More generally, there are three classes of solutions for this model whose 4D metric is 
maximally symmetric [2^, and all of them have vanishing vacuum energy. The solution 
displayed above is the simplest one, which corresponds to putting branes of equal tension at 
the antipodes of the the compact space. A second class of solutions has branes of unequal 
tension, in which case the metric no longer factorizes, but instead is warped. The third 
class of solutions, which are also warped, is more exotic: it has nonconical singularities at 
the antipodal points whose stress energy does not correspond to simple brane sources. In 
particular, the Trr and Teg stress tensor components are nonvanishing for these backgrounds. 



1. The argument against self-tuning 



In ref. 41[, it was argued that, rather than dis play ing self-tuning, the SUGRA model is a 
clear example of Weinberg's no-go theorem. Ref. 41[ dimensionally reduces the model with 
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two equal-tension branes to 4D gravity coupled to the dilaton and the radion ijj, which 
enters the metric ansatz in the form 

_ ^^f. ^ M^^e-^^ {dr^ + sin^ r de^) (10) 

The potential which results for the two scalar fields is found to have the form 

Viij, 0) = Mg-^ e'^^ (J^ e-'"^ - 2Ml e"'^^ + 2A6^ (11) 

where a2 = — (f), ai = 2ip + (p, and a is related to the deficit angle induced by the branes, 
a = 1 — r/(27rM|). This potential looks different from the dilaton potential in because 
in the latter, the Einstein equations have already been imposed, which fixes ^p in terms of 
and Ag, whereas this has not been done in (lllj). Ref. ^| goes on to observe that, starting 
from some configuration in which V{il!, 0) = 0, any change in the brane tensions, hence in 
the value of a, will make Vltp, 0) nonzero, and induce a runaway potential for the field a2- 
If the 4D potential (fTTj) really captures the essential features of the model, then clearly 
there is no self-tuning; however this ansatz describes only the unwarped solutions. If self- 
tuning occurs, it might be through a transition from an unwarped to a warped solution, 
especially an exotic nonconically singular one, which cannot be described by the ansatz 

(uni). 

More generally, the proposal of involves large extra dimensions. 



with Mq^ ~ 0.1 mm, where the Kaluza-Klein mass gap is of the same order of magnitude 
as the mass of the stable direction ai in the potential (fTTjl . The argument of j4l| has thus 
been called into question, on the basis that there is no justification for ignoring the role of 
light KK modes while keeping other modes that are of the same energy in the effective 4D 
description. This objection is similar in spirit to the previous one, in that mild warping of 
the background metric can be thought of as a coherent state of KK modes. 

2. The argument for self-tuning 

One reason for suspecting a self-tuning mechanism in the Salam-Sezgin model is the 
observation by that the only solutions with maximal symmetry in the 4D part of the 
metric are those with vanishing curvature in 4D, that is, Minkowski space. These solutions 
fall into two classes, one of which has only conical singularities (or possibly no singularities). 
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It has the warped hne element and dilaton profile 



^ f2 + f2 \^-^ [1 + f^/flY (1 + f7ff)2 j 

^0 + In fS^) (12) 



0. 



where 



■•"—^ '•-IF (13) 

In the special case where Tq = (z.e., Age''^'' = (3'^e~'^° /2), we recover the unwarped 
spherical football-shaped solutions, which have branes of equal tension at the antipodal 
points. For Tq 7^ r^, the tensions of the branes are related to each other by 

(1 - At^/M^) (1 - At./M^) = rl/fl = (14) 

where is an integer, related to the quantization of gauge field flux on the compact manifold. 
The most general solutions are more complicated; they are given by 



0(^) = (po + 4:lnW + 2X37] (15) 



with 



,4 _ ^2 /e<^o/3 coshAi(r7 - t]i) 



= — 



a 



Ai V 2A6 cosh A2(r7 - 772) 



and A3 = yAl — Af. When A3 = 0, this reduces to the special solutions above. 

One can see the nonconical nature of the singularities when A3 7^ by performing a 
coordinate redefinition [4^. Writing dr = aW^drj, in the limit r] ±00 

dsl = dr^ + r^''^d9^ (17) 

where 

_ ±2A3 + A2 + 3Ai 

±2A3 + 5A2-Ai' ^ 

As soon as A3 7^ 0, one finds that a± 7^ 1 and the metric near the branes no longer 
corresponds to that of a simple conical singularity. Rather, we find a geometry more like 
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that of a trumpet, where the circumference of closed circles around the singular points is 
not linearly proportional to their radius. It is the existence of such solutions which has 



41|. Indeed it is not excluded that under 



given hope of circumventing the argument of 
a perturbation to the branes' tensions, the system would relax to a static solution with 
nonzero A3 rather than move to a runaway solution with purely conical singularities. It is 
the presence of such behaviour that we will be looking to confirm or negate in the analysis 
that follows. 



C. The need for thick branes 



The solutions we have summarized aboved only concern pure tension branes and are 
therefore not suited for the analysis we wish to carry out, where the energy density on the 
branes is time dependent, and its equation of state must be different from p = —p. In order 
to study this more general case, it is necessary to regularize the singularity at the brane 
positions by giving the branes finite thickness. The reasons for this were explained in detail 
in We will only recall here the main line of reasoning. 

For a metric of the form 

ds^ = -n\f,t)dt^ + a\f,t)dx^ + f{f,t){df^ + f^de^) (19) 

CO dimension-two branes will appear in the stress energy tensor as (p,p) x 

6{f)/[2TTf'^{f,t)f]. There must be terms in the Einstein tensor that provide the necessary 
delta functions to match the ones from the stress-energy tensor. Since the two dimensional 
delta function is given by 

6{f) = 2nfV^ ln(f) (20) 

the delta function terms in the Einstein tensor will come from terms like ln{n), ln(a) 
and V^ln(/). In order to accommodate general equations of state, the first two of these 
must be nonzero near the branes, which means that the warp factor must scale like f"*^*) as 
f — i> 0. As was done in to treat the Einstein-Hilbert case, we will be smoothing out 
this mildly singular behaviour by giving the brane finite thickness, allowing us to find well- 
behaved solutions with arbitrary equations of state, and thus arbitrary time dependence of 
the stress-energy tensor, on the branes. 
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III. THICK CODIMENSION-TWO BRANEWORLD IN 6D SUPERGRAVITY 

For the reasons that have just been exposed, we need to regularize the branes by giving 
them nonzero thickness in order to consistently study their cosmology in situations where 
their stress-energy is not simply constant. We will be following the formalism of fl^, where 
such a construction was used to conclusively rule out the presence of self-tuning in the 
context of six dimensional Einstein-Hilbert gravity. The only new ingredient here will be 
the presence of a dilaton field. 

We start with the most general dynamical metric having axial symmetry in the extra 
dimensions 

f^s2 ^ ^)2^^2 ^ ^(^^^ ^^)2^^2 ^ ^^^^ ^Ydr'^ + c(r, tfde'^ + 2E{r, t)drdt (21) 

and perturb its components in the following manner 

n(r,t) = e^°('^)+^^ a(r, t) = ao(t)e"^"('^)+^i('^'*); 6(r, t) = 6o(t)e^''^''^+^'^'''*^; 

c(r,t) = Co (t)e^«(")+^^ E{r,t) = Ei{r,t). (22) 

The zeroth order solutions are for simplicity taken to be ones corresponding to zero-tension 
branes, and the perturbations represent the effect of small time-dependent brane stress- 
energies. We also expand the dilaton and gauge potential as 

e'^^'-'*) = </.o(t)e^°('-)+^^('"'*^ Ag{r,t) = Af\r) + (23) 

The stress energy tensor is given as 

n = tl + sl + s^l (24) 

where represents contributions from all bulk fields, and and s^^l represent the brane 
content. Since we are assuming that the branes are only present at the perturbative level, 
their stress energy can be written as 

s\ = -pe{ro - r) - S{ro - r).Fo(t); = pe{ro - r) - S{ro - r)J^^{t)- 

si = plB{ro-r)] Sq = pgOiro - r); s[ = p[ 0(ro - r); s* = ^(ro - r); 

sjf- = -p^6{r - r^) - 5{r - r^)J^^{t); s^l = p^6{r - r^) - 5{r - r^)J^^{t); 

sj^ = pj^6{r - r^); sjg = p/gd{r - r^); = p^l 6{r - r^); sj^ = pj;.6{r - n) 

(25) 
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where 6{r) are Heaviside step functions. Our branes are therefore represented by cores of 
finite radii tq and located around the poles of an axially symmetric compact internal 
space. For simplicity, we will assume and tq are such that both branes have the same 
thickness at the level of the unperturbed background geometry. 

Notice the inclusion of one-dimensional delta function terms in the brane stress-energy. 
In H, 

these came from expanding the step function in a Taylor series, treating the time 
dependence of the thickness as a perturbation. Here, since the background tension is assumed 
to vanish, such terms do not appear unless we put them in by hand. The end result however 
can be shown to be exactly the same: the one- dimensional delta functions effectively encode 

n n 

the time dependence of the brane thickness. (See and j4a] for discussions on this point). 
These functions are not assumed to be the same at both branes (i.e. we will not be demanding 
that JFq = JF^), so that the above assumption that the background thicknesses are equal is 
not overly constraining. 

Finally, we will assume that time derivatives of 0{p) perturbations are of 0{p^^'^), which 
is implied by the usual law for conservation of energy p ~ (d/a)p ~ p^^"^. 

A. Background solutions 

With the above ansatz for the metric and matter content, working in coordinates where 
Bo{r) = one can show that 

= _^g-4Ao(r)+Co(r)+<Ao(r-)_ ^26) 

dr 

Since at this level there is no stress-energy on the brane to induce warping or singularities, 
the only possible solution is the perfectly spherical static "football" solution jj^l, where 

Ao{r) = (27) 

Mr) = 00 (28) 
e^o(0 = (29) 

(Note that the coordinate r defined here differs from the one in section IIA by the rescaling 
r — i> kr. In these coordinates, 6 ranges over (0, 27r) even in the presence of a deficit angle, 
while r of course ranges over (0, Tc/k).) Substituting these into the background equations of 
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motion, one finds 

^o{t)v{%) = ^ (30) 

2v9o(t)co(t)^Oo(i)^ 

— -| = ^o{t)v{^^) + ^ (31) 

One can readily check that these imply that 



3^-^ (32) 



and that the time dependent functions 6o(^))Co(t) and (po(t) obey 



6o(t) = Co(t) (33) 

y,o(t) = , ^ . (34) 

/2t;($o)co(t)2 



Notice here a key difference with the Einstein-Hilbert case, where there is no dilaton. In 
this case, we would have (po{t)e'^°v{^o) Aq and P"^ /[{po{t)e'^°] — > and we would instead 
find that &o(^) and co(t) must be constants. This difference can be understood as follows. In 
the Einstein-Hilbert case, since the radion is stabilized we find in our perturbative analysis 
that the scale factor of the internal space is static. In the SUSY case, there are two degrees 
of freedom: the radion and the dilaton. Only one linear combination of the two is stable, 
while the other corresponds to a fiat direction. It is this mode which appears here as a new 
dynamical entity. 

B. Perturbed equations of motion 

Before writing down the equations of motion to linear order in the perturbations, it is 
useful as was done in to find a set of variables which are invariant under the gauge 
transformation 

t^/(r,t); r^g{r,t). (35) 

One convenient set is 

Z = N[~A[- W = 3A[ + N[- X = ^~B,-%C,- 

P ~ '^*ti P ~ Pt ~ ~^ ^*t^ Pr ~ ~^ ~^ ~^ ^*t^ 

p5 = sl + p6 = Sg- si + s4 - s^, (36) 
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in terms of which the perturbed equations of motion become 



W - C',W = co(t)2^ 



Z' + c^z 



Oo 

ao 



Y' - C'qY - e'^°/?e^''(") {W - <P\) 
C'^W , ^/2^;(<fo)e-^«M 
Co(t)^ 



ao Co 



Co 



+ 



Co 

Co 



+ 



(P + P) 







+ 



coityMi 

ao 



2t;(<l>o)e'^« 



Co(t)2M| 



+ 



+ ^ + 2^^ 
Co Vco 



1 

M, 



p', - C',p, = 

U' -Y - pe^^e^'^^'^X = 

-coitm 



Pi 



ao 4 M| Co 

■p + 3^(p + + ^(2p + 2p5 + pe) = ft + « 
ao Co 



Oo 

ao 



+ 



1 Co _ Oo Co ^ p5 

3 Co ao Co 3Mg 



(37) 

(38) 
(39) 

(40) 
(41) 
(42) 
(43) 

^'i) (44) 
(45) 

(46) 



C^X' 2Xp^2v{^o)e'^' ^ 3C^W 



2v{<l>o)e 



-Co(r) 



coity 



1 



2co(t)^ 



-Y 



coityMi 

ao V ao / Co c^o Cq 3 V Co / 



(47) 



As was done in we assume all the p's and p's are functions of time only. Also, in order 
to simplify calculations and because in the absence of a specific microscopic model for the 
brane matter content we are free to choose pQ as we wish, we will assume that 



P6 



V, t) = e{ro - r)e2^°WP6(t) + 0{r - (vr/fc - ro))e2^°Wp,6(t). 



(48) 



C. Solutions to the perturbed equations of motion 

The preceding system of equations may seem daunting, but is actually quite straight- 
forward to solve. One simply uses eqs. ^7\i andpij) to solve for Z{r,t),W{r,t) and 
P5(r, t) respectively. Once these solutions are found, one can use eq. (pUj) to solve for 0i(r, t) 
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and then either eq. (jlUj) or eq. (jH^ to solve for Y{r, t). Then X(r, t) can be solved using eq. 
(jlTj), while eq. (jl^j) solves for p[(r, t) and either eq. (jl^ or eq. pijl solve for U{r,t). 

One must then impose appropriate boundary conditions, namely that all solutions are 
regular at the poles, and smooth across the core/bulk boundary except where indicated 
otherwise by the presence of the one-dimensional delta function terms in the stress-energy 
tensor. 

In the interest of brevity we will only present here the general solutions for W{r,t) and 
0i(r, t) in the bulk, which will be required for our argument in the following section. 



0f"'*^)(r,t) = ^ J^r^'it) + J^r^'it) ln(l - cos(A;r)) + J^^^'it) + irnzlZ^v— ^ (t) 



1 



(bulk),,. . cos{kr) (bulk) 



+ 



2k 

Coit)'^ ( f Co\ ^Hq ^ ^CqCIo n, f 1 (^{bulk),,^ ^(bulk) . 



2k 



- _ 3- + 3- - - 3 p - -- ij^r >it) + rr 'it) 

k \Vco/ Oo CqQq \aoJ I 2k ^ 



ln(sin(/cr)) 
(50) 



where J^i{t) are constants of integration set by boundary conditions. The interested reader 
will find solutions to the other gauge invariant variables in the appendix. 



IV. PRESENCE OF THE A3 / SOLUTIONS 

As we have explained above, the hope is that starting from a static A3 = solution 
the system will naturally evolve to a static A3 7^ solution when the brane tensions are 
perturbed. In order for our results to confirm or rule out this idea, we must verify that our 
perturbative solutions actually include the A3 7^ case. We will now show this explicitly. 

We work in the convenient gauge where -Bi(r, t) = 0, which means that effectively gf^r = 1- 
In this gauge we can write the metric as 

ds^ = -n{r)dx^dx^, + dr"^ + c{rfde'^ (51) 

and it was shown in that 



^3 = ^c(r) 



2n(r)V(r)' + Anirfnir)'] . (52) 



In our perturbative language, the background solution has a static dilaton and no warping 
so that A3 vanishes identically at this order. If we wish to prove that our perturbative results 
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include the nonconical solutions with A3 7^ 0, we have to show that at 0{p), the combination 
on the right hand side of eq. does not vanish. For static perturbations, we can write 
eq. ^ to C(p) as 

5A3 = ^e^°M[20i(r)' + iy(r)]. (53) 

Substituting our bulk solutions into this expression and assuming staticity (i.e. ao(t) = 
1, co(t) = 1, p = —p and p=K = — p*), we find that 

1 



= Ulkli ' cos(fcro))' iP.it) - P.6(t)) . (54) 



This result tells us that the solutions with nonconical singularities at the poles are related 
to the presence on the branes of stress-energy along the extra dimensions (recall that we 
defined — sj; = 9{ro — r)p^{r, t) = 9{ro — r)e^'^°^'^^VG{t) and sjg — sj^ = 9(r — r^)p^r, t) = 
6{r — r^)e'^^°^^'^V^Q{t) ). This result is compatible with the discussions of 43| and |45[. 

Having shown that our solutions include the case of interest, we can now address the 
question of whether the model can solve the cosmological constant problem through self- 
tuning. 



V. COSMOLOGY AND THE QUESTION OF SELF-TUNING 

A. Boundary conditions and the Friedmann equations 

In braneworld models, one derives the Friedmann equations on the brane(s) by impos- 
ing appropriate boundary conditions. In the model we are studying, these imply that all 
functions are regular at the origin, and smooth across the brane/bulk boundaries. Since 
all the gauge invariant quantities we have defined only appear as first r derivatives in the 
equations of motion, the functions themselves must match, but not their first derivatives. 
The only exception is 0i(r, t), for which both the function itself and its first derivative must 
be matched across the boundaries at r = tq and r = n/k — tq. One must also take into 
account the one-dimensional delta functions in the brane stress-energy tensor, which affect 
the following junction conditions 

limX(ro - e,t) - X(ro + e,t) = ^-^^t^co{tf:Fo{t) (55) 
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limX(r^, + e, t) — X(r^, — e, t) 



tan(fcro) ^ ^^^2 



Co(t)'^.(t) 



p[(ro,t) = -J-o(t) - 2^^o(t) 

Co 

Co 



(56) 
(57) 
(58) 



Solving for all variables and boundary conditions, one is led to the Friedmann equations, 
and conservation of energy 

^ (1 — cos(A;ro)) , , (1 — cos(A;ro))'^ 



ooV 



8Mi 



-{p + p + p* +p*) 



72A;2Mi 



*6) 



3 Co Vco/ 
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Consistency between these four equations implies that 
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(63) 



Note that the constant of integration Qi which appears in eq. ()63|) can be shown to be 
proportional to the perturbation to the flux. 



B. Effective four dimensional quantities 

As was discussed in the quantities we have used so far are not the ones a four 
dimensional observer on the brane would identify as the energy density and pressure. Indeed, 
to define the effective four dimensional quantities, we must integrate the 6D quantities over 
the thickness of the brane, 

S^^^l = 27r / dr 6(r, t) c(r, t) S^^^^l (64) 
Jo 
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which with our choice of gauge {grr = 1) leads to 

pro 

Jo 

= 27r / co(t)'e^°('-VW^^r - 27rco(t)'e^»(^°) JoW (66) 
Jo 

Jn/k—ro 
Jn/k—ro 

We also define the 4d Newton constant as 

= Af| = 2,r jf'" c„(()^e°»MAf|<ir (69) 

We come now to the main result of this paper, the Friedmann equations on general 
CO dimension-two branes in six- dimensional supergravity. Writing these in terms of the ef- 
fective four dimensional quantities defined above, we find 
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(75) 

We see a major difference with the Einstein-Hilbert case studied in here the 4d New- 
ton constant is time-dependent at first order in the perturbations, whereas in the Einstein- 
Hilbert case such corrections would only appear at O(p^). Although we will not focus on this 
point for the rest of the paper, we could use results of fifth force experiments and measures 
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of the time variation of Newton's constant to strongly constrain the dynamics of the variable 
Co(t). In fact, one can show from the equations presented above that in the "Pg = 'P*6 = 
limit, the equations can be derived from the effective action 

1 



S = d xy— ^ 



IGttG 



1 

2C^ 



+ Cm (76) 



by making the identifications Co(t) = -\/C, V"(C) = — 327r^ ^vnikr^GQxjk and G jC^it) = G^^t) 
In other words, the model in that case corresponds to Brans-Dicke gravity with u = 1/2 
and a flat potential, a theory which is clearly ruled out by experiment. Note that this is 
consistent with the arguments of |41|, since this is precisely the limit where non-conical 
singularities are absent, and corresponds to the case that was studied in that work. 



When one does include non-zero Vq and P^g terms, we see from eqs. ()7H|74|1 that con- 
servation of energy is violated from a four dimensional perspective, and it is thus hard to 
imagine how one might write down an effective four dimensional action. It should be noted 
that this should not come as a surprise, since it is precisely this type of intrinsically extra di- 
mensional effect which motivated this approach as a plausible way to circumvent Weinberg's 
no-go theorem and arrive at a self-tuning solution to the cosmological constant problem. 



C. Absence of self-tuning 

Using eqs. (f7T|) - (f73jl . we can now categorically answer in the negative the question of 
whether or not the model self-tunes to a static solution under perturbations to the brane 
tensions. 

Indeed, it was already known that in the absence of the solutions with nonconical sin- 
gularities, one needs to fine-tune the brane tensions with the flux from the bulk gauge field 



in order to obtain static solutions 



S Q, 0, Q, Q 



The hope, as explained above, was 



that including the solutions with nonconical singularities would allow the model to naturally 
evolve from a static solution with conical singularities to a static solution with nonconical 
singularities when one perturbs the tensions. We see clearly now that this cannot happen, 
for two reasons. 

First, nonconical singularities are only present when the functions Vq and "P^g are nonzero. 
However, these have no dynamical equation of motion; they are free parameters in this 
model. Therefore, even if it were possible to use them in order to cancel a change in the 
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brane tensions, this would have to be done by hand and would thus represent an arbitrary 
fine-tuning of the model's parameters. This would be made all the more unnatural by the 
fact that it would have to be a time-dependent fine-tuning if it were to hold across multiple 
phase transitions on either brane. 

Second, it is not even possible to use the nonconical singularities to cancel perturbations 
to the brane tensions and get new static solutions. This can be seen by setting p^'^^ = — p*^^) 
and pi^^ = —pi^^ in the above equations, and setting all time derivatives to zero. We then 
find that eq. (fTTj) necessitates a relation between p^'^\ p^f' and Qi that, when plugged into 
eq. (|7H|l tells us that = —7^*6- In other words, one must tune the brane tensions in order 
to obtain a static solution regardless of whether we have conical or nonconical singularities. 
Moreover, if we do have nonconical singularities, we must perform an additional fine-tuning 
of the functions Vq and 'P^.e sourcing them in order to keep the solutions static. Therefore, 
we have conclusively shown that no self-tuning mechanism for the cosmological constant is 
present in the type of six dimensional supergravity we have considered here. 



We should note that although this final result agrees with the conclusions of j4l|, the 
analysis we have performed here was far from superfluous, since we have shown that self- 
tuning fails not only for the simple case considered in that work, but in the much more 
general case which includes warping and nonconical singularities. Indeed, the solutions we 
find are qualitatively different from the ones in when Vq and "P^g are nonzero. One 
cannot simply derive them from an effective four dimensional action since they correspond 
to a theory where the ususal law of conservation of energy does not hold (see eqs. ()7Hf74|l ). 

One final point of interest is the fact that one can see from our results that we cannot set 
only one of and Cq constant without imposing unnatural constraints on the brane stress- 
energy tensor components (see eq. ()75|) in the case Cq = est for example). They must either 
both be constant (which requires the fine-tuning mentioned in the previous paragraphs), or 
both be time dependent. This explains why previous solutions always singled out solutions 
with Minkowski 4-space over de Sitter or anti-de Sitter space: the initial assumption was 
always made that the internal space was static! Indeed, rather than saying that once the 
radion and dilaton are stabilized, solutions with a static external space are required by the 
equations of motion, a more correct statement would be that choosing the large dimensions 
to be Minkowski is actually part of the tuning one has to do in order to obtain a static 
radion and dilaton. 
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D. Dropping axial symmetry: a loophole? 

Although our results appear quite conclusive regarding the absence of a self-tuning mech- 
anism for the cosmological constant in six dimensional supergravity, one should be keep in 
mind which assumptions were made in deriving this result, to see if relaxing one of them 
might in fact lead to a different conclusion. 

The solutions we have found are more general than those that had so far been presented 
in the literature because they allow for general equations of state for the brane content and 
they include the possibility of warping and of nonconical singularities in a fully dynamical 
context. There is however one restriction we have imposed on our metric ansatz which could 
be relaxed: axial symmetry of the internal space. 

In a recent paper explicit solutions have been found with an arbitrary number of 
3-branes on a compact two-dimensional internal space by dropping the requirement of axial 
symmetry. The analysis in [46| concludes that the only solution with two branes is one with 
equal tension branes located at the poles of a spherical internal space, but this result is 
derived by excluding the possibility of warping and is restricted to Einstein-Hilbert gravity. 
It is interesting to ask what might happen in the supersymmetric case if we were to add 
the possibility of breaking axial symmetry in response to a change in the brane tensions, 
in addition to the possibility of nonconical singularities and warping. While there is no 
compelling reason to believe that allowing the branes to shift their positions would lead to 
self-tuning, we feel it is nonetheless an avenue worth investigating in future work, since it 
could conceivably lead to interesting and perhaps unexpected effects. 

VI. CONCLUSION AND OUTLOOK 

In this paper, we have studied the cosmology on a co dimension-two brane in six- 
dimensional supergravity. We achieved this by solving the dynamical field equations lin- 
earized around a static background. In order to deal with the mildly singular behaviour 
of the warp factor at the position of the branes when the equation of state of their matter 
content is different from that of pure tension, we regularized the branes by giving them a 
nonzero thickness. 

Our results show that there is no self-tuning of the cosmological constant to zero in such 
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a setup. Rather, any change to the branes' tensions leads to expansion in both the external 
space and the internal space. In fact, we find that the internal and external space must 
always either both be static (which necessitates fine-tuning) or both be time dependent, 
but that we cannot have one be static while the other evolves without imposing unnatural 
relations between the brane stress-energy tensor components. This explains why previous 
solutions, which always assumed a static internal space from the start, singled out static 
solutions in the large dimensions. 

Our solutions include both nonconical singularities and warping of the metric, and we 
emphasize that this represents a new result, since no dynamical treatment of such a model 
including both these possibilities had been carried out to date. The only assumption we 
made is that of axial symmetry of the internal space, and although it would be interesting 
to look for solutions where this assumption is dropped, we see no reason to believe that it 
would lead to self-tuning. 

Although we do not believe that 6D supergravity theories offer a solution to the cosmo- 
logical constant problem, such models with large extra dimensions can nevertheless provide 
a novel source of dynamical dark energy (whose size is naturally of the right magnitude once 
the cosmological constant is set to zero) as well as accelerator signals which would be visible 
at the LHC j4J] , and they therefore merit more study along these lines. 
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APPENDIX A: GENERAL SOLUTIONS TO THE PERTURBED EQUATIONS 

OF MOTION 

We present here the general solutions to the perturbed equations of motion (jH7j) - (jTfj) . We 
only present the solutions in the core located in the region < r < tq. Bulk solutions have 
the same form, without the terms coming from the brane stress-energy, while solutions in 
the other core have the same form with the terms from the brane stress-energy replaced by 
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the corresponding terms for this second brane (p* instead of p for example). The solutions 
for the integration "constants" J-i{t) will not be given explicitly, but can be found in each 
region of the model through the imposition of boundary conditions, as explained in the main 
text. 
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We have not written down U{r, t) explicitly since it can be obtained trivially from the 
algebraic equation (j44|) once all other functions are given. 
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